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Abstract. - Phase diagram and pattern formation in two-dimensional Ising model with coupling 
between order parameter and lattice vibrations is investigated by Monte-Carlo simulations. It 
is shown that if the coupling is strong enough (or phonons are soft enough) a short-range order 
exists in disordered phase for a broader temperature interval. Different types of this short-range 
order (stripe- like, checkboard-like, etc.) depending on the temperature and model parameters are 
investigated. With further increase of the coupling, a reconstruction of the ground state happens 
and new ordered phases appear at low enough temperatures. 



Introduction. — Traditional views assume that un- 
der the thermodynamic equiUbrium conditions a system 
should either be homogeneous or consist of macroscopi- 
cally large domains of homogeneous phases. It appears 
frequently, however, that equilibrium or very long-lived 
metastable states occur with a nanoscale or mesoscale het- 
erogeneity (for a general review of these phenomena in var- 
ious physical and chemical systems see, e.g., Refs. [1-4]). 
It is commonly accepted now that the formation of the 
mesoscale heterogeneity is a result of frustrations in the 
system which can result from either geometric factors [4-7] 
or competing interactions, the long-ranged forces such as 
Coulomb or dipole-dipole interactions being of primary 
importance [4,8-14]. 

"Heterophase" fluctuations in metallic alloys [15] pro- 
vide us interesting examples of intrinsic nanoscale inho- 
mogeneities. In such cases, a system behaves like an en- 
semble of nanosize particles of one phase embedded into 
the matrix, e.g., so-called athermic w-phase in bcc host 
which is observed in some Ti- and Zr-based alloys [16], as 
well as in Cri_a;Ala; [17]. This peculiar structural state 
leads to strong anomalies of electronic properties [18] and 
ultrasound attenuation [19]. The "athermic w-phase" is a 
term to describe rather some short-range order in atomic 
positions than real phase since it never exists in the bulk 
(atomic positions in this "phase" are intermediate between 
bcc structure and true w-phase existing in Ti and Zr un- 



der pressure). It would be not surprising to observe strong 
short-range order in the close vicinity of a second-order 
phase transition. On the contrary, these heterophase fluc- 
tuations sometimes arise in a broad temperature and con- 
centration domain. Nature of this state is still unknown 
(see discussions in Refs. [7,15]). 

A well-pronounced short-range order or nanoscale inho- 
mogeneity in a broad temperature interval are often ob- 
served in magnetic alloys with a strong coupling between 
magnetism and lattice (or chemical composition), such as 
Cu-Mn alloys [20,21] or Fe-Ni Invar alloys [22,23]. As 
mentioned above it is a tendency now to connect intrinsic 
inhomogeneities in various systems with long-range inter- 
actions. One may expect therefore that a long-range char- 
acter of elastic deformations is solids could be relevant for 
the problem under discussion. 

To clarify a role of elastic-mediated interactions in pos- 
sible pattern formation we have investigated the simplest 
model of Ising order parameter coupled with phonons at 
the square lattice. Basing on the results of Monte Carlo 
simulations for this model we show that, indeed, this cou- 
pling can result in a formation of various nanosize-scale 
structures. 

Model and computational details. The two- 
dimensional Ising model [24, 25] is a prototype, exactly 
solvable model of order-disorder phase transitions in mag- 
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netic systems, ordering alloys, etc. To be specific we 
will use further terms "magnetic" and "spins" to describe 
the ordering phenomena under consideration. To describe 
the effects of magnetoelastic (spin-lattice) interactions, we 
proceed with the Hamiltonian 



(1) 



where Aj^ = 1 if i,j are the nearest neighbors and is equal 
to zero, otherwise, Uj are atomic displacement vectors and 
spin variables Si = ±1. The spin-lattice coupling is taken 
into account via coordinate dependence of the exchange 
parameters Jij = J(Ri + Uj — Rj — Uj), Rj being square 
lattice vectors. Assuming that the displacements are small 
the exchange parameters can be written in the linear ap- 
proximation 



u 



(2) 



where riij is the unit vector in direction of Rj — Hj. 

Substituting Eq.(2) into Eq.(l) and assuming periodic 
boundary conditions we obtain 



(3) 



where ^"J^ = KnfjU^jAij . The dependence of the ex- 
change parameters on interatomic distances results in the 
last term in right-hand-side of Eq.(3) 



-^T — '^''^■ijJij^i^j 



(4) 



describing forces acting on atoms due to spin redistribu- 
tion. These forces initiate the displacements of atoms into 
new equilibrium positions R -|- uq determined by the ex- 
pression 

<(R0 = EG"MR-'^)^/ (5) 
j 

where R^ = Rj — TLj, Gc«/3(Ry) is the static lattice 
Green's function determined by the equation 



(6) 



Replacing the variables u = uq -|- g one can represent the 
partition function Z as 



Si 



(7) 



where the effective spin Hamiltonian 

Heff = J E 4SiSj + J E P^"Ga|s{niJ)P^ (8) 



and /3 = 1/T is the inverse temperature. 

Due to harmonic approximation for the potential en- 
ergy of atomic displacements and linear approximation for 
the magnetoelastic coupHng the phonon and spin subsys- 
tems turn out to be totally separated after the change 
of variables. The statistical properties of the system is 
completely determined by the Hamiltonian H^ff which 
depends on spin variables only. The last term in the right- 
hand-side of Eq.(8) reads 



lj2J2''?k'^ikSiSkGc.0{nij)n''^iJ'jiSjSi 



(9) 



ij kl 



and describes indirect long-range spin-spin interactions via 
lattice distortions. 

To investigate equilibrium properties of the spin system 
with the Hamiltonian (8) we have carried out the Monte 
Carlo simulations for the square lattice with 100 x 100 
sites and periodic boundary conditions, using a standard 
Metropolis algorithm [26] . To provide the Gibbs distribu- 
tion P ~ exp(— /3ile//) we have performed up to 10* fiips 
for each lattice spin. 

We are interested in the two-dimensional spin system 
embedded into three-dimensional elastic medium. For 
two-dimensional continuum the Green's function is patho- 
logical, with logarithmic growth at large distances. In- 
stead, we use the expression 



a/3 



Sap Si{2TTR/a) 



A-KIJL 



R 



(10) 



valid in the framework of quasi-continuum approach [27] 
with isotropic Debye model for the phonon spectra. Here 
Si{x) is the integral sine function [28], a is the lattice 
parameter, /i is the shear modulus which can be ex- 
pressed [27] in terms of the force constants The 
Green's function (10) has asymptotic behavior GapiR) oc 
1/R at large distances. To speed up computations, 
we use its truncation by multiplying expression (10) by 
exp(-(i?/L)*). The length L and ph onon-induced inter- 
action cut-off radius were chosen at ■5tli and 13th neigh- 
bors, respectively; these values ensure the convergency and 
stability of computational results. 

For the bare exchange parameters the nearest- 
neighbor (nn) interactions have been taken into account. 
The choice of nn parameter | | = 1 determines the energy 
units. 

Results of the Monte Carlo simulations. — In the 

square-lattice Ising model with nn interactions the phase 
transition takes place at T = Tc = 2.264. One can see 
(Fig. 1) that the true ground state is reached for approx- 
imately 10^ spin fiips per spin whereas for shorter Monte 
Carlo runs metastable configurations such as domain walls 
appear. Further we will use by default this number of fiips. 
Probably it is not enough for a close vicinity of Tc due to 
critical slowing down but we will not discuss this region. 

The effects of long-range interactions due to magnetoe- 
lastic coupHng are crucially dependent on dimensionless 
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Fig. 1: The pattern of spin distribution at T = 1 (a.b) and T = 3.5 (c.d) without spin-lattice coupling for the AFM case. 
Simulation results with 10"^ (a,c) and 10^ (b,d) flips per spin are presented. Black and white regions indicate the spin-up and 
spin-down, respectively. 




Fig. 2: The pattern of spin distribution at T = 1.5 (a,b,c) and T 
at g = -10.0 (a), Q = -4.5 (a), Q = -2.5 (a), Q = -10.0 (a). 



3.5 (d) with the spin-lattice coupling for FM ground state 



coupling constant Q — {J')^ / fijf . Nontrivial patterns 
in a broad enough temperature interval both below and 
above Tc arise at \Q\ > 3.5; typical spin configurations are 
shown in Fig. 2. In this case we have found a long-range 
or short-range order of checkboard or stripe types with 
staggered spin-up and spin-down regions. These patterns 
appear when MC simulations start from both random and 
from regular initial spin distribution and, therefore, are 
not metastable but, rather, equilibrium states. Pictures 
of the long-range order are quite similar for the case of 
FM (Jf = -1) and AFM (Jf = 1) Ising model. 

Since these new types of ordering result from the cou- 
pling between magnetic and elastic subsystems it is natu- 
rally to expect its manifestations also in the distribution 
of the displacement field (Fig. 3) . For weak magnetoelas- 
tic coupHng the lattice distortions are created by domain 
boundaries in metastable configurations and spread over 
large distances due to a slow decay of G{r). Fig. 3a shows 
the corresponding distribution calculated by the pertur- 
bation theory in the lowest order in Q. Spin patterns at 
large Q are inevitably related with the lattice distortion 
chessboard patterns visible in Fig. 3b. Contrary to pre- 
vious case, deformations are mutually cancelled at large 
spatial scale. 

To understand character of new types of ordering let 
us consider energies of different local spin configurations 





Fig. 3: Distribution of atomic displacements iti 
case at T = 1 without spin lattice coupling {Q = 0) (a) and at 
T = 3.5 with the spin-lattice coupling (Q = —9.0) (b). Black 
regions correspond to Ux+Uy < and white ones to Ux+Uy > 0. 
The results are presented for 10^ flips per spin. 



(Fig. 4). For ferromagnetic ordering (Fig. 4a) the force 
(4) P" = and therefore the magnetoelastic contribu- 
tion to the total energy vanishes. On the other hand, 
the configuration Fig. 4b corresponds to the largest local 
value of P" and, thus, to the maximal magnetoelastic en- 
ergy gain. This configuration is a minimal structural block 
of stripe and checkboard structures. Direct total energy 
calculations for various spin configurations (sec the Ta- 
ble) demonstrate that at Q ^ —3.5 simple FM or AFM 



p-3 



I. K. Razumov et al. 



Table 1: Energies (per atom) of different spin configurations; 
FM, AFM, means the FM and AFM nearest-neighbour inter- 
action; "c" rnxn labels checkboard ordering (see, e.g. Fig. 2b) 
with elementary white and black rnxn rectangles; "s" mxn la- 
bels diogonal stripes with elementary steps rn and n in a; and 
y directions; "R" labels random spin configuration. 





FM 


AFM 


Q 


-4 


-3 


4 3 


FM 


-4.0 


-4.0 


4.0 4.0 


AFM 


4.0 


4.0 


-4.0 -4.0 


R 


-1.8 


-1.0 


-1.8 -1.0 


c2x2 


-5.3 


-3.0 


-5.3 -3.0 


c3x3 


-4.2 


-2.9 


-1.6 -0.3 


c2x3 


-4.8 


-3.0 


-3.5 -1.7 


s2xl 


-4.7 


-2.6 


-4.7 -2.6 


s3xl 


-4.2 


-2.9 


-1.5 -0.3 
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Fig. 4: Local spin configurations forming different patterns. 






Fig. 5: A phase diagram for ferromagnetic case. 1, 2, and 3 
labels FM phase, checkboard ordering, and disordered (para- 
magnetic) phase, respectively. 



structure become energetically unfavorable and the new 
ordered ground state occurs. In this case the checkboard 
ordering is preferable and the stripe structure has slightly 
higher energy. 

With the temperature increase, the long-range order 
is destroyed. Wc estimate critical temperatures of these 
transitions calculating the i?-dependence of pair correla- 



tion functions of the corresponding order parameters. The 
results are presented in Figure 5. With the growth of Q 
the temperature of transition from FM (AFM) and para- 
magnetic phase decreases. This is not surprizing since 
the magnetoelastic coupling destabihzes these spin states. 
With further increase of Q there is a transition between 
check-board (or stripe) phase and the paramagnetic state, 
the transition temperature growing with Q. 
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Fig. 6: Space dependence of the correlation function < SoSx > 
in X (or Y) direction for Q = (solid curves 1, 2, 3, and 4) and 
Q = —2.0 (dashed curves 1', 2', 3', and 5). The temperatures 
are equal to 0.9 Tc (curves 1, 1'), 1.1 Tc (curves 2, 2'), 1.2 Tc 
(curves 3, 3'), 1.3 Tc (curve 4), and 1.65 Tc (curve 5). 

Figure 6 shows the evolution of the spin correlation func- 
tions for IQI < 3.5 (that is, at the transition between 
regions 1 and 3 at the phase diagram 5) with the tem- 
perature increase. One can see that the magnetoelastic 
coupling leads to stronger short-range order above Tc (c.f. 
curves 2 and 2', 3 and 3', respectively), which survives up 
to a relatively higher temperature, in comparison with the 
original Ising model. In particular, the short-range order 
for Q = —2.0 and T = 1.65Tc (curve 5) is stronger than 
for Q = and T = 1.3Tc (curve 4). 

The short-range order in the region 3 near the triple 
point tends to formation of stripe-like structures (Fig. 
2b, c), although this configuration is not the most ener- 
getically preferable. The appearance of this structure is 
connected with its higher entropy in comparison with the 
checkboard structure. 

Discussion and conclusions. — Our computational 
results are summarized schematically in the phase diagram 
shown in Fig. 5. For small enough coupling constant 
IQI < 3.5 the behavior of the system is qualitatively sim- 
ilar to that of the standard Ising model. In this case the 
main effect of the magnetoelastic coupling is suppression 
of Tc due to destabilization of FM or AFM states (Fig. 6). 
Usually this happens if Tc is strongly suppressed by soft 
fiuctuations, such as in quasi-two-dimensional Heisenberg 
magnets [29] or frustrated three-dimensional Heisenberg 
magnets with competing interaction [30] . In those cases Tc 
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is much smaller than a typical energy of nearest-neighbor 
interaction J whereas the short-range order in a broad 
temperature interval survives up to the temperatures of 
order of J. We illustrate here an alternative mechanism 
of occurrence of the short-range order due to coupling with 
a second subsystem providing long-range effective interac- 
tions. 

Thus, in this regime the magnctoclastic interactions are 
favourable for the short-range order in paramagnetic phase 
but do not result yet in the nontrivial pattern formation. 
The latter happens at Q > 3.5 where exotic spin config- 
urations (checkboard or stripe states) become favourable. 
Especially, near the triple point at the phase diagram (Fig. 
5) there arc many different states with approximately the 
same free energy and therefore a formation of complicated 
inhomogeneous states could be expected. At further tem- 
perature increase we reach a conventional paramagnetic 
(random) state without short-range order. This scenario 
reminds the behavior observed at the "melting" of stripe 
domains in magnetic films [11], our "pattern" regime being 
similar to "tetragonal phase" in the latter case. The main 
formal difference between these two problems is two-spin 
character of long-range dipolc-dipole interactions versus 
four-spin character of long-range phonon-mediated inter- 
action considered here. 

Despite a simphcity of our model, it demonstrates a 
rather general feature which we believe can be relevant for 
discussions of real systems demonstrating a strong short- 
range order above Tc in a broad temperature interval. 
The effect arises at large enough coupling constant \Q\ 
which means either unusually strong interaction between 
the subsystems (large J') or soft phonons (small shear 
modulus). At least, in some case, such as Cu-Mn alloys the 
shear modulus tends to zero at some critical composition 
or temperature [21] so \Q\ can be, in principle, arbitrarily 
large. It would be interesting to analyse the problem of 
heterophase fluctuations in alloys mentioned in the Intro- 
duction from this point of view. 
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